Anisotropic spin relaxation in graphene 
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Spin relaxation in graphene is investigated in electrical graphene spin valve devices in 
the non-local geometry. Ferromagnetic electrodes with in-plane magnetizations inject 
spins parallel to the graphene layer. They are subject to Hanle spin precession under 
a magnetic field B applied perpendicular to the graphene layer. Fields above 1.5 T 
force the magnetization direction of the ferromagnetic contacts to align to the field, 
allowing injection of spins perpendicular to the graphene plane. A comparison of the 
spin signals at B = and B = 2 T shows a 20 % decrease in spin relaxation time for 
spins perpendicular to the graphene layer compared to spins parallel to the layer. We 
analyze the results in terms of the different strengths of the spin orbit effective fields 
in the in-plane and out-of-plane directions. 

PACS numbers: 72.25.Rb,72.25.Hg 



The discovery of the anomalous quantum HaU 
effect in graphene [l|, Q triggered an avalanche of 
theoretical and experimental work on this new sys- 
tem. Spintronics is one of the fields which hasgreat 
expectations for this material. Spin qubits [3| and 
many other spintronic devices based on graphene 
could become available due to the fact that in in- 
trinsic graphene spins are expected to relax very 
slowly [j, [1, • The reason behind this is the 
low hyperfine interaction of the spins with the car- 
bon nuclei (only 1 % of the nuclei are C13 and have 
spin) and the weak spin-orbit (SO) interaction due 
to the low atomic number. 

Recent experiments show spin trans por t in 
graphene up to room temperature [1,0, [13, 11. 12|. 
with spin relaxation lengths of 2 /im and relaxation 
times around 150 ps Such relatively short re- 
laxation times suggest an important role of SO in- 
teraction. There are two relevant mechanisms for 
SO interaction in graphene [l3|. In the EUiott- 
Yafet (EY) mechanism, spin scattering is induced 
by electron (momentum) scattering from impuri- 
ties, boundaries and phonons. The Dyakonov- 
Perel (DP) mechanism results from SO terms in 
the Hamiltonian of the clean material. Here elec- 
trons feel an effective magnetic field, which changes 
in direction every time the electron scatters to 
a different momentum state, resulting in random 
spin precession. In principle the two mechanisms 
can be distinguished by their different dep endence 
on the momentum scattering time r |13| . In our 
experiments in graphene we are not able to change 
r considerably, making the distinction between 
Elliott- Yafet and Dyakonov-Perel mechanisms dif- 
ficult. However we can obtain valuable informa- 



tion about the SO interaction by investigating the 
anisotropy of spin relaxation. First we note that 
the transverse {T2) and longitudinal (Ti) spin re- 
laxation times are expected _to_ be the same for the 
parameters of our system 



13| 



Therefore, as in 
metals, a single spin relaxation time T = Ti — T2 
can be used. However, due to the 2-dimensionality, 
T can have a different value for injected spins par- 
allel (Ty) or perpendicular {T±) to the graphene 
plane. For example, if the SO interaction is of 
the Rashba or Dresselhaus type then the SO ef- 
fective fields are exclusively in the graphene plane 
and calculations show that this should result in 
anisotropic spin relaxation in which T±_ — 1/2 Ty 
[l3|. On the other hand, if the SO effective fields 
pointing out of the graphene plane dominate, we 
expect Ti >> T||. Here we will directly compare 
the spin relaxation times in the parallel and per- 
pendicular directions, measured under identical ex- 
perimental conditions. 

Our experiments are performed using the four 
terminal "non-local" technique (Fig. la). Here 
the charge current path can be fully separated 
from the voltage detection circuit. The non-local 
technique is less sensitive to device resistance fluc- 
tuations and magnetoresistanccs (such as Hall ef- 
fects), as compared to the standard two-terminal 
spin valve technique. This allows the detection of 
small spin signals, in our case as small as a few 
ToSl (Fig. 2a). Fabrication of the devices is done 
as in Ref. 8. Using the "Scotch tape" technique 
[ist graphene layers were deposited on an oxidized 
(500 nm) heavily doped Si wafer. Calibrations by 
Raman spectroscopy in combination with optical 
microscopy and atomic force microscopy show that 
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FIG. 1: Spin transport in graphene a) A SEM pic- 
ture of a single layer of graphene contacted by 6 cobalt 
electrodes (sample A). Spins traveling a distance of 5 
/im, from cobalt electrode F2 to F3, are probed us- 
ing the 'non-local' geometry. The voltage circuit {F3- 
graphene-F6) is completely separated from the current 
circuit (Fl-graphene-F2), b) Hanle type spin preces- 
sion experiment, the magnetization of the spin injector 
F2 is set antiparallel to the magnetization of spin de- 
tector F3. Spins are injected parallel to the graphene 
plane, c) Application of a strong external magnetic 
field (~1.4T) perpendicular to the graphene layer re- 
sults to injector and detector magnetizations aligned 
parallel to the external magnetic field. Spins are in- 
jected perpendicular to the graphene plane, d) Spin 
precession in case of parallel (||, black curve) and an- 
tiparallel (II, grey curve) magnetizations. 



our samples are single graphene layers. We evap- 
orate a thin layer of aluminum (6A) on top of the 
graphene layer at 77 K and let it oxidize using 
pure O2, to form an AI2O3 barrier. These bar- 
riers very likely contain pinholes nevertheless 
spin injection efficiencies of 10 % have been ob- 
served. Conventional electron beam lithography 
and e-beam evaporation of 50 nm of Co (at 10~^ 
mbar) are used to define the ferromagnetic cobalt 
electrodes. The electrodes have different widths to 



assure different switching fields 1J|. The experi- 
ments are performed at a temperature of 4.2K and 
we use magnetic fields up to 4.5 Tesla. A standard 
a.c. lock-in technique is used with currents in the 
range 1-20 fiA. 

Spin precession measurements are performed on 
two samples (graphene width W = 1.2 fim) for 
electrode spacings L = 5 ^m (sample A), 0.5, 2 
and 4/im (sample B). To perform the Hanle spin 
precession experiments we first apply a magnetic 
field in the y-direction to prepare the magnetiza- 
tions of the electrodes in a parallel or antiparallel 
orientation (Fig. la). Then this field is removed 
and a i?-field in the z-dircction is scanned (Fig. 
lb) [1^ . An example of the resulting spin pre- 
cession is depicted in Fig. Id (sample B), for the 
parallel and antiparallel magnetizations of the spin 
injector and spin detector cobalt electrodes. The 
spins are injected parallel to the graphene plane 
and are precessing while diffusing towards the spin 
detector situated at a distance i = 4 /im from the 
injector. At ^ 0.2 T the average precession 
angle is 180 degrees, resulting in a sign reversal 
of the spin signal. The magnitude of the signal 
at Bz = OT (0.2 fl) is small compared to the sig- 
nals measured in our previous work, which was in 
the order of 2 for these spacings 0]. This is 
directly related to the low contact resistances Rc 
(1-2 kO) found between the cobalt electrodes and 
the graphcne/Al203, which are a factor 5 to 10 
smaller than in Ref. @. In this study, the con- 
tact resistance Rc is equal or smaller to the typical 
square resistance of the graphene layer R^q and this 
results in the reduction of the injection/detection 
efficiencies and also provides an extra path for spin 
relaxation at the ferromagnetic contacts [3| ■ This 
is taken into account later in the fitting of the spin 
precession measurements with solutions of the 1- 
dimensional Bloch equations [l3| which describe 
the combined effect of diffusion, precession and 
spin relaxation in the system. From the fit, the 
intrinsic spin relaxation time in graphene can be 
extracted. The relative importance of the contacts 
is given by the parameter R = W ■ Rc/Rsq, where 
W is the width of the graphene layer. If the spin 
relaxation length is in the fim range then the model 
shows that for R >> 10^^ m the contacts do not 
induce extra spin relaxation. On the other hand 
for R << 10~^m the amplitude (A) of the spin 
signal has quadratic dependence in R (A ^ R^) ( 
M. Popinciuc et al., in preparation). 

In our experiments a background arises in the 
non-local resistance which has a quadratic depen- 
dency on Bz due to orbital magnetoresistancc ef- 
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FIG. 2: Anisotropic spin relaxation at high electron 
density n a) n = 2.1 10^^ cm~^. Initially, spins are in- 
jected parallel to the graphene plane having the mag- 
netization of the spin injector set parallel (gray line) or 
antiparallel (black line) to the detector (at a distance 
1/ = 5 ^m). From the fits (dashed line) we extract the 
diffusion constant D and the relaxation time Ty. A 
magnetic field of 1.4 T is needed to align the magne- 
tization of the cobalt electrodes out of their easy mag- 
netization axis, in this case the spins are injected per- 
pendicular to the graphene layer having a spin relax- 
ation time T± . The spin relaxation T± is 19 % smaller 
to Ty . The small decrease in the nonlocal signal found 
between 1.5 and 2 T is attributed to a background due 
to orbital magnetoresistance effects. The same exper- 
iment has been performed on sample B for b) L = 4 



fim [n = 3.5 10^^ cm ^), c) L — 2 fim [n 
cm"^) and d) L = 500 nm (n = 3.5 10^^ cm 
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fects, in particular around the Dirac point. At high 
electron densities this background is usually small. 
We therefore start by applying a gate voltage on 
samples A and B to allow us to investigate the 
spin dynamics at a high electron density n~ 3.0 
10^^m~^(Fig. 2). From the fitting procedure both 
the diffusion constant D and the spin relaxation 
time T|| can be obtained. For L = 2, 4 (Sample 
B) and 5/im (Sample A) we obtain a diffusion con- 
stant of 3- 10~^TO^s~^ and spin relaxation times of 
Tij = 60 ps for sample A up to 90 ps for sample B, 



corresponding to spin relaxation lengths A=^£T|| 
of 1.4 up to 1.8 /im, respectively, comparable to the 
values found in Ref. [8| . Increasing the magnitude 
of Bz results in a rotation of the magnetization of 
the cobalt electrodes out of the plane, towards the 
magnetic field direction. This can already be seen 
in Fig. Id where the rotation of the magnetization 
induces an asymmetry of the spin signal at 0.5 T. 
A magnetic field of 1.4 to 1.8 T is needed to fully 
align the magnetization of the cobalt electrodes in 
the z-direction (Fig. Ic) [l4, 17|. This is a special 
situation as injected spins are now perpendicular 
to the graphene layer and will relax with a time 
Tj_ which is not necessarily the same as T|| . If the 
anisotropy in the spin relaxation is large then the 
amplitude of the non-local spin signal at B = T 
should be very different from the signal at ^^1.8 
T. The high sensitivity of this method comes from 
the fact that the amplitude of the spin signal de- 
pends exponentially on A. For example, for L = 
5 /xm and A = 1.4 /xm a decrease by 10 % in A 
due to anisotropic spin relaxation should result in 
a 40 % decrease in spin signal amplitude. In this 
example a 10 % decrease in A corresponds to a 20 
% (100% • (Til - Ti_)/Ti_) decrease in T. In Fig. 
2 the decrease in the magnitude of the spin sig- 
nal for L = 2, 4 and 5 /xm corresponds to a spin 
relaxation time T± being 20 % smaller than T||. 
Clearly, our devices show anisotropic spin relax- 
ation in graphene at high electron densities. Of in- 
terest is to investigate if the same conclusion holds 
for spins injected in graphene at the charge neu- 
trality point. For L = 2, 4 and 5 /xm, close to the 
Dirac point, orbital magnetoresistance effects in- 
duce a large background, increasing quadratically 
in Bz- This background is not only monotonic in- 
creasing, it also contains non-periodic fluctuations 
function of Bz with an amplitude equal or 
larger than the spin signal. This effect, in combi- 
nation with the large suppression of the spin signal 
amplitude at the Dirac point, for i = 2, 4 and 5 
/um, does not allow us to investigate in precision 
the spin anisotropy. However, we were able to per- 
form the experiment for the L = 0.5 /im spacing in 
which the spin signal is relatively large. (Fig. 3c) 
Here, application of a gate voltage of -76.5 V (Fig. 
3a) allows us to investigate the spin dynamics at 
the Dirac point. Clearly, the non-local resistance 
at 2T is smaller than the resistance at Bz = OT 
showing similar anisotropic spin relaxation behav- 
ior as for high electron densities. 

We can now estimate the effective magnetic field 
which the electrons feel, assuming the Dyakonov- 
Perel mechanism. The electron scatters to a dif- 
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FIG. 3: Anisotropic spin relaxation at the Dirac point 
for L = 500 nm a) Gate voltage dependence of the 
graphene resistance between electrodes F2 and F3 (Fig. 
lb). The charge neutrality point is found at Vg « -75 
V. b) The nonlocal resistance for spins injected parallel 
to the graphene layer (black line). This is defined as the 
difference in the signal obtained for injector and detec- 
tor magnetizations set parallel (Ratl.tt) ^^'^ the signal 
for injector and detector set to antiparallel (Rjv_t,||). 
Our model (gray line) which takes into account the 
finite contact resistance gives a qualitatively good fit 
to the data for Vj<30V. c) Application of a magnetic 
field B perpendicular to the graphene plane gives sim- 
ilar results as found in Fig. 2. A spin relaxation for 
spins injected perpendicular to the graphene layer is 
found to be smaller than the spin relaxation of the 
spins injected parallel to the layer. 



ferent momentum state after a time r which re- 
sults in a precession angle of the spin Aw = WpT. 
Here, cUp is the precession frequency of the spin. 
The number of scattering events necessary to in- 
duce an angle of 27r is y^T/r. Using T = 100 ps 
and T ~ 30 fs we obtain ujp ~ lO^^s^^. Therefore, 
in the Dyakonov-Perel mechanism, the precession 
frequency corresponds to a magnetic field of about 
5 Tesla. 

We now check if the transverse (T2) and lon- 
gitudinal (Ti) relaxation times for spins in the 
graphene plane are the same. For this, A {=^/DTi) 
is extracted from the spin signal dependence L 
(Fig. 4). Care has to be given to the strong 
suppression of the spin signal amplitude as we ap- 
proach the Dirac point, found at 4.2 K and as well 
at room temperature. Our model takes into ac- 



count the spin relaxation at the contacts and fits 
qualitatively well the data (Fig. 3b). Earlier work 
did not show this strong effect due to the fact that 
in those samples the contact resistances were large 
enough to not to influence the spin dynamics. In 
Fig. 4 we present the amplitude of the spin signal 
of sample B as function of the electrode spacing L 
(= 0.5 , 2 and 4 /xm) for different values of R. In 
the same figure we present the length dependence 
in the signal expected from our model. At high 
electron densities n (Fig. 4a) the model gives a 
spin relaxation of 1.5 ±0.2/im. This value is com- 
parable to the values found from spin precession 
measurements (Fig. 2) performed at similar val- 
ues of n (and i?), proving that Ti ~ T2. The effect 
on A when we approach the Dirac point is stronger: 
our model gives A = 1 /xm for n ~ 2.0 lO^^cm"^ 
(Fig 4 b) and A = 0.8 /xm for n ~ 1.0 lO^^cm'^ 
(Fig. 4c). 

Summarizing, we use a non-local measurement 
geometry to investigate anisotropic spin relaxation 
in graphene. Although our experimental accuracy 
do not allow us to confirm that Ti and T2 are 
equal, they allow us to measure with accuracy the 
difference between T±_ and T||. At high electron 
densities, a decrease up to ~40 % in the spin sig- 
nal is found for injection of spins perpendicular to 
the graphene layer compared to injected spins par- 
allel to the graphene. This corresponds to a spin 
relaxation Tj_ almost 20 % smaller to T||. This 
spin anisotropy is expected for a 2-D system were 
spin-orbit fields in plane dominate the spin relax- 
ation. As a next step we suggest to investigate 
the dependence of Tl and T\\ on momentum scat- 
tering T to establish the relative role of the EUiott- 
Yafet mechanism compared to the Dyakonov-Perel 
mechanism. 
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FIG. 4: Dependence of the spin signal as function of 
electrode spacing L for spins injected parallel to the 
graphene plane, a) The gate voltage Vg is set in such 
a way that for L = 0.5 , 2 and 4 we have the same 
value of i? = 2.2 10"*^ (n ~ 5.0 10^^ cm-^ D = 0.04 
m'^s"^). We obtain a relaxation length of A = 1.5 ± 
0.2 which is similar to the length extracted from 
the spin precession measurements (see Fig. 2). Model 
calculations for A = 0.7^im and A — 5/im are shown for 
comparison, b) Moving towards the Dirac point using 
a gate voltage in such a way that we decrease the R 
value to 1.0 10"^ (n ~ 2 10^^ cm-^ D = 0.03 m^s"^) 
has a strong influence in A, as it decreases to 1 and 
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